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We present the discovery of a new gravitational lens system with two compact 
radio images separated by 0.701±0.001 arcsec. The lens system was discovered in the 
Cosmic Lens All Sky Survey (CLASS) as a flat spectrum radio source. Both radio 
components show structure in a VLBA 8.4 GHz radio image. No further extended 
structure is seen in either the VLA, MERLIN or VLBA images. Hubble Space Telescope 
(HST) WFPC2 images in F555W and F814W show two extended objects close to the 
radio components, which we identify as two lens galaxies. Their colours and mass- 
OO ! to-light ratios seem to favour two late-type spiral galaxies at relatively high redshifts 

Q\ . (zd <; 0.5). Faint emission is also detected at positions corresponding to the radio 

images. 

A two-lens mass model can explain the observed VLBA structure. The best fit 
model has a reduced x 2 °f LI. The relative positions of the VLBA subcomponents 
are reproduced within 0.08 mas, the flux density ratios within 0.19. We also reproduce 
the position angle and separation of the two VLBA subcomponents in A and B within 
the observational errors, which we consider strong evidence for the validity of the lens 
model. Moreover, we find a surface density axis ratio of 0.74±° j° for the primary 
lens (Gl), consistent with the surface brightness axis ratio of 0.69 ± 0.15. Also, the 
surface density position angle of (4. 91^2 4)° °f Gl compares well with the (—6 ± 13)° 
position angle of the surface brightness distribution. The errors indicate the 99 per 
cent confidence interval. 
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1 INTRODUCTION 

In the last few years gravitational lensing has proved useful 
not only in the determination of cosmological parameters, 
such as the Hubble constant (e.g. Refsdal 1964, 1966) and 
the cosmological constant (e.g. Kochanek 1996), but also in 
the study of the mass distribution in the universe and the 
mass distribution of lensing galaxies. To obtain a sample 
of gravitational lens systems, relatively unbiased compared 
to optical lens surveys (Kochanek 1991), which suffer from 
seeing effects and dust obscuration, two large radio surveys, 
the Jodrell-Bank VLA Astrometric Survey (JVAS; Patnaik 
et al. 1992; King et al., in preparation; Wilkinson et al., in 
preparation) and the Cosmic Lens All Sky Survey (CLASS; 
Myers et al., in preparation), were set up. Together these 



surveys targeted ~12,000 flat spectrum radio sources with 
flux densities larger than 25 and 200 mjy for CLASS and 
JVAS, respectively. All sources were observed with the Very 
Large Array (VLA) in A-array at 8.4 GHz with 0.2 arsec 
resolution. Objects that showed signs of multiple compact 
components, or structure that could be due to lensing, were 
listed for further high resolution radio observations with 
MERLIN. Those objects still exhibiting compact structure 
in the MERLIN image were subsequently observed with the 
VLBA to confirm their identification as a lens system and 
sometimes HST to observe the optical emission of the lens 
galaxy and lens images. 

In the following sections we give a detailed description 
of B1127+385, a newly discovered gravitational lens system. 
In Section 2 we describe the radio observations. In Sections 
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Table 1. VLA 8.4 GHz (1), MERLIN 5 GHz (2), MERLIN 1.7 
GHz (3) and VLBA 8.4 GHz (4) astrometry and flux densities 
for B1127+385. The integrated flux densities from the VLBA 
observations are given. The VLBA positions of Al and Bl are 
ll h 30 m 0.099 s , +38°12'3.091" and ll h 30 m 0.157 s , +38°12'3.232" 
(J2000), respectively. The errors on the flux densities are ~10 per 
cent. All errors are 1<j. 
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3 the optical HST observations are presented. In Section 4 
we present a lens model based on the image positions and 
flux density ratios from the VLBA observations. In Section 
5 we summarise our results and conclusions. 



2 RADIO OBSERVATIONS 

2.1 VLA and MERLIN observations 

B1127+385 was observed on 1995 August 14 with the VLA 
in A-array at 8.4 GHz as one of the ~10,000 flat spectrum 
CLASS sources. The image shows two compact components 
separated by ~0.7 arcsec (Fig. 1). The 92-cm (0.327 MHz) 
WENSS (de Bruyn et al., in preparation) flux density is 
12 ± 3 mjy, establishing a slightly inverted radio spectrum. 
This implies that both components, with similar flux den- 
sity at 8.4 GHz, most likely have a flat or inverted spectrum. 
This immediately made B1127±385 a strong lens candidate, 
because a change alignment within 0.7 arcsec of two unre- 
lated compact flat spectrum radio sources is less than 10~ 6 . 
Thus there is only a probability <1 per cent of finding such 
a change alignment in the sample of ~10,000 flat spectrum 
sources. No indication of variability has been found sofar. 
A MERLIN 1.7 GHz long-track observation on 1996 Febru- 
ary 13 showed two compact components. Subsequently, a 
snapshot observation with MERLIN was made at 5 GHz in 
1996 December. These observations show only the two com- 
pact components A and B (Fig. 1). The VLA and MER- 
LIN observations were reduced in AIPS and mapped by the 
Caltech package DIFMAP (Pearson et al. 1994; Shepherd 
1997). The flux densities and positions were determined in 
DIFMAP by fitting Gaussian models to both components 
simultaneously. The results are listed in Table 1. The spec- 
tral indices of the components are —0.08 (A) and 0.05 (B) 
{S„ cx v a ) respectively between the VLA 8.4 GHz and MER- 
LIN 5 GHz flux densities, and -0.05 (A) and -0.11 (B) 
respectively between the VLA 8.4 GHz and MERLIN 1.7 
GHz flux densities. The errors on these spectral indices are 
~ 0.15. The compactness and similarity in spectral index 
of both radio components underline the lens candidacy of 
B1127+385. But high resolution radio observations and op- 
tical follow-up are necessary to secure its lensing nature. 



Table 2. The separations and position angles (north to east) of 
components A2 (B2) with respect to components Al (Bl). The 
flux densities are those of (Al, A2) and (Bl, B2), respectively. 
The errors on the flux densities are ~ 10 per cent. All errors are 
la. 

Sep. (mas) P. A. (°) 5 8 . 4G Hz (mjy) 
A 2.20 ±0.15 16 ±4 10.5,3.2 
B 0.97 ±0.15 155 ±9 7.9,2.9 

2.2 VLBA observations 

VLBA 8.4 GHz observations were made on 1996 Novem- 
ber 4. Phase referencing was used, switching between 
B1127±385 (4 min integration) and the strong nearby JVAS 
phase-reference source B1128±385 (Patnaik et al. 1992; 2 
min integration) for a period of 3.5 hours. The map of 
B1128±385 shows a sub-mas unresolved point source at 8.4 
GHz. Fringe-fitting (Thompson, Moran & Swenson 1986) 
was therefore performed on B1128+385 and the solutions 
were directly transferred to B1127±385, which is ~11 ar- 
cmin west of B1128±385. All data reduction was performed 
in AIPS and mapping was done in DIFMAP. The data was 
uniformly weighted. The resulting map resolution is 1.7 x 0.7 
mas (PA of 41°). 

Components A and B show clear evidence for the pres- 
ence of substructure (Fig. 1). Model fitting in DIFMAP 
shows that two axisymmetric Gaussian components can well 
represent the substructure in images A and B. The posi- 
tions of the Gaussian components and their flux densities 
are listed in Table 2 for their best fit (minimum \ 2 )- The 
separation between Al (A2) and Bl (B2) is 700.7 (700.0) 
mas and the position angle of the line from Al to Bl is 
78.5°. The flux density ratios between Al and Bl, and A2 
and B2 are 1.3 and 1.1, respectively. The integrated flux 
density ratio of A over B is 1.3, consistent with the 8.4 GHz 
VLA (1.3) and 5 GHz MERLIN (1.3) flux density ratios. 
The large position angle difference of ~ 139° between the 
subcomponents in A and B is expected if A and B are lensed 
and given opposite parities (Schneider, Ehlers & Falco 1992). 



3 HST OBSERVATIONS 

Hubble Space Telescope (HST) exposures of B1127+385 in 
the filters F555W (V) and F814W (J) were taken on 1996 
June 21, using the Wide Field Planetary Camera (WFPC2). 
The exposures were taken on the PC chip (45.5 mas pixel -1 ) 
and the exposure times in V and I band were 700s±300s and 
2 x 500s, respectively. A standard reduction was performed 
on both images. The /-band image is shown in Fig. 2. The I- 
band image shows two clear emission peaks within ~1 arcsec 
distance from the radio components. Except for a bright 
nearby galaxy 8 arcsec south, no other galaxies are seen near 
B1127±385. Because the absolute astrometry of the HST is 
poorly matched (offsets of ~1 arcsec) to the more accurate 
VLBI astrometry, we cannot 'blindly' overlay the optical and 
radio maps. However, the contour plot of the optical /-band 
emission, convolved to 0.1 arcsec (Fig. 3), clearly shows two 
bright (Gl and G2) and one fainter emission feature. If we 
assume that the radio component A is associated with the 
optical emission (~ 9a peak) west of Gl, we find that also 
radio component B is associated with an emission feature 
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Figure 1. Upper left: VLA 8.4 GHz snapshot of B1127+385 taken 1995 August 14. The image has a resolution of 0.2 X 0.2 arcsec 
and shows two compact components (B is east, A is west) separated by 0.7 arsec. Contours arc at (—3, 3, 6, 12, 24) X 0.367 mjy per 
beam. Upper right: MERLIN 5 GHz observation taken 1996 December. The map has a resolution of 0.06 X 0.04 arcsec (PA of 44°). 
Contours are at (-3,3,6,12,24,48,96) X 0.146 mjy per beam. Lower left: VLBA 8.4 GHz image of B1127+385 component A, taken 
1996 November 4. Contours are at (—4, 4, 8, 16, 32, 64) X 0.090 mjy per beam. Lower right: Idem of component B. Contours are at 
(-4.5,4.5,9, 18,36,72) X 0.078 mjy per beam. The VLBA maps have a resolution of 1.7 X 0.7 mas (PA of 41°). 



(~ 5cr peak). Although there appears to be a slight offset 
between the optical emission and the radio position of B, this 
could be due to the poorer signal-to-noise or the extended 
nature of the optical emission near B. 

Both Gl and G2 are extended, suggesting that both 



are galaxies. Photometry and relative astrometry were per- 
formed on Gl and G2 in I band, and only photometry in V 
band (Table 3). The V — I colour indices of Gl and G2 are 
1.9 and 2.0 mag, respectively. The separation is 0.60 arcsec 
and the position angle of the line G1-G2 is 120°. 
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Figure 2. HST /-band image of B1127+385. The two radio com- 
ponents are associated with two emission features west (A) and 
cast (B) of Gl (see also Fig. 3). North is up, east is left. The area 
shown is 3.1 arcscc X 2.6 arcsec. 



Table 3. HST relative astrometry (7 band) and photometry 
for B1127+385. Component Gl is located at ll h 30 m 0. 1726 s , 
+38°12'1.903" (according to the STSDAS/metric routine in 
IRAF). All errors are la. 

Att(mas) A<5(mas) I (magn.) V (magn.) 
Gl 0±5 0±5 22.5 ±0.1 24.4 ±0.2 

G2 521 ± 5 -299 ± 5 23.5 ± 0.1 25.5 ± 0.2 








£^ 

Figure 3. Contour plot of B1127+385. The HST image has been 
convolved to 0.1 arcsec to bring out the optical emission more 
clearly. The two markers show the radio positions of A and B, 
where A has been placed on the centre of the optical emission as 
determined by a Gaussian fitting procedure. The contours indi- 
cate (3, 6, 9, 12, 15, 18, 21) X rms noise in the image. 



4 MODELLING 

In this section we present a model that reproduces the ob- 
served properties of B1127±385. We use a Singular Isother- 
mal Ellipsoid (SIE) mass distribution (Kormann, Schneider 
& Bartelmann 1994) to describe the lens galaxies. We as- 
sume a smooth FRW universe. If not mentioned otherwise, 
all errors indicate 99 per cent confidence intervals. 

From the VLBA observations we obtain 10 constraints 
(8 from the image positions and 2 from the flux density 
ratios). The two source positions give 4 free parameters and 
the mass model gives 3 (velocity dispersion, surface density 
(SD) axis ratio and position angle). The number of degrees 
of freedom (NDF) is therefore 3. 



4.1 Single lens mass model 

Initially we try a model consisting of a single SIE galaxy. We 
place the mass distribution on the surface brightness (SB) 
centre of Gl, determined by fitting a 2D-Gaussian profile to 
it. The position of Gl relative to the optical emission feature 
associated with radio component A is (—0.228", —0.038"), 
with a la error of ~5 mas in both x and y. The parities of 
components A and B are taken as -1 and ±1, respectively. 
Choosing ±1 and -1 for A and B, and letting the centre 
of the SD distribution move freely, give in all cases (both 
for the single and double lens case) unsatisfactory models 
(X 2 £ 5). 

Using the image positions and flux density ratios of the 
radio components from Tables 1 and 2, we project rays back 
on the source plane through the lens. For the two image 
pairs - A1&B1 and A2&B2 - we simultaneously minimize 
the distance between the back-projected rays and the differ- 
ence between observed and model flux density ratios (Kayser 
1990). We allow for a la error of 0.1 mas in the relative x and 
y distances between the positions of the two back-projected 
rays and their average position. A la error of 0.15 is allowed 
for the flux density ratios. When the model has converged 
sufficiently, we use the average source positions to calculate 
the image positions in the lens plane. These are subsequently 
used to calculate a \ 2 from the mismatch with the observed 
image positions. The resulting mass model parameters for 
minimum \ 2 (lens plane) are listed in Table 4 (model I) . 



4.2 Double lens mass model 

Although Gl is the primary lens, G2 cannot be neglected as 
it lies close to Gl and is only ~1 mag fainter. We therefore 
extend the model by placing a Singular Isothermal Sphere 
(SIS) at the position of G2 — (-0.749", -0.337") relative 
to component A, with a la error of ~5 mas in x and y. 
We choose a SIS for G2, because G2 is only an external 
perturber of the primary lens Gl and as long as there is no 
need to complicate the model (poor x 2 ) one should keep the 
mass model as simple as possible. Using a SIE for G2 would 
add two extra free parameter (axis ratio and position angle) 
and make the mass model less constrained. 

The velocity dispersion of G2 is fixed at 0.79 • cP 1 , us- 
ing L oc <T|| (e.g. Faber & Jackson 1976) in combination 
with the ~1 mag difference between Gl and G2 and assum- 
ing they have similar mass-to-light ratios. We repeat the 
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X 2 minimization procedure. The resulting model parame- 
ters are listed in Table 4 (model II). Model II reproduces 
the relative positions of all four images to within 0.08 mas, 
compared to 0.16 mas for model I. And although the num- 
ber of degrees of freedom does not change between models 
I and II, the reduced \ of model II is significantly smaller 
than that of model I. A reduced \ 2 > 4-5 corresponds to a 
probability of 4 • 10 -3 and model I can therefore be rejected 
as an appropriate model with 99.6 per cent confidence. So 
adding G2 improves the mass model significantly. A faint 
image will be formed to the south-east of G2, which can be 
removed by a very small core radius (0.01") for G2, without 
changing the model parameters at any signifcant level. The 
critical and caustic structure of model II is shown in Fig. 4, 
where G2 has been given a core radius of 0.01". 

If Gl and G2 are spiral galaxies, the relation between lu- 
minosity and velocity dispersion is L oc a 2 ' 5 (TuUy & Fisher 
1977), hence of 2 « 0.69 • of 1 . The best model then has 
a slightly increased \ 2 of 1.6, still a considerable impove- 
ment over model I. The resulting model parameters deviate 
by less than 3 per cent in velocity dispersion and axis ra- 
tio from model II. The position angle of Gl becomes —5.3°, 
even in better agreement with the observed surface bright- 
ness position angle (see below). However, all values are well 
within the errors determined by Monte-Carlo simulations for 
model II (see Section 4.4 and Table 4). 

Using a SIE for G2 would have the most influence on 
the position angle and axis ratio of the mass model of Gl. 
However, the close agreement between these parameters and 
their observed values (see next paragraph), as well as the 
small resulting minimum x 2 , indicate that adding extra free 
parameters is not necessary. 
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Figure 4. The critical curves (thick lines) and caustics (subpanel) 
of model II in Table 4. The star (subpanel) indicates the model 
position of the source as seen in the source plane. The solid con- 
tours depict the images at the positions of the radio components, 
if the source surface brightness is gaussian with a FWHM of 5 
mas. The lens galaxy Gl is indicated by the cross-haired cllips, 
G2 by the cross-haired circle. The dashed contours indicate the 
time-delay surface for = 0.5, z s = 1.5, Ho=50 km s — 1 Mpc -1 , 
Qm = 1 and Q\ = 0. 



4.3 Surface density versus surface brightness 

The SD axis ratio of Gl, (b/a)^ 1 = 0.74 (model II), is only 
slightly larger than the SB axis ratio of 0.69 ± 0.15 that we 
find from fitting a 2D-Gaussian profile to the SB distribution 
of Gl. The same Gaussian fit gives a position angle of (— 6± 
13)°, close to the SD position angle of 4.9°. Although the SB 
profile of Gl is most likely not Gaussian, the position angle 
and axis ratio inferred from a 2D-Gaussian fit will give a 
good indication of the value for these parameters. 

Strong evidence that B1127+385 is a gravitational lens 
system is given by the expected centre of the SD distribution 
of Gl. In Fig. 5 the 99 per cent confidence contour of the 
central SD position of Gl is plotted. The two circles indicate 
the radio components A and B. If we assume that the two 
faint optical emission features are associated with A and B, 
we find that the optical emission peak of Gl (cross) falls 
perfectly inside the 99 per cent confidence contour. In other 
words, the position of the SB centre of Gl relative to the 
faint emission feature west of it coincides with the position 
of the SD centre of Gl relative to radio component A. This 
suggests Gl and G2 are indeed two lens galaxies and the 
faint optical emission features (Fig. 3) are associated with 
the radio components. 



4.4 Monte-Carlo simulations 

To investigate the overall reliability of model II, we per- 
formed Monte Carlo simulations. We minimize y 2 for 10,000 
models, where we add Gaussian distributed errors (la) to 
the relative image positions (0.1 mas), flux density ratios 
(0.15), galaxy positions (5 mas) and velocity dispersion ra- 
tio between Gl and G2 (0.10). 

4-4-1 Mass model parameters 

The parameter probability density distributions of all mod- 
els with x 2 < H-3 (99 per cent confidence interval) are 
shown in Fig. 6. The figure also shows the 99 per cent con- 
fidence interval (shaded) of the observed SB axis ratio and 
position angle. We see that the probability distributions of 
both the SD axis ratio and position angle have considerable 
overlap with these shaded regions. Hence, the axis ratio and 
position angle of the luminous matter agree well with those 
inferred from the SD distribution of Gl. The mean values of 
the recovered position angle and separation of Al relative 
to A2 are (18 ± 7)° and 2.0 ± 0.2 mas and for Bl relative to 
B2 (163 ± 3)° and 1.1 ± 0.2 mas, where the errors are the 
rms values of the parameter probability distributions. These 
recovered model parameters compare well with the observed 
values listed in Table 2. We consider this as strong evidence 
for the validity of the lens model. 

To investigate the ratio (cr|| >G2 /cr|| > Gi): we calculate the 
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minimum \ f° r a range of this ratio, as shown in Fig. 6. 
The two horizontal lines indicate the 90 per cent and 99 per 
cent confidence intervals of the \ 2 distribution. The shaded 
region indicates the 99 per cent confidence interval for the 
ratio (<T||,g2/o"|| i gi)) determined above. The dot gives the ra- 
tio we determined from the ~1 mag luminosity difference. 
We see that this ratio lies well below the 99 per cent con- 
fidence level, which shows that the ratio expected from the 
mass model agrees with that determined from the F-J rela- 
tion. Also the ratio from the T-F relation agrees well. 

4-4-2 Time delay 

The predicted time delay between components A and B 
is 1.51±g;|| ■ /a, s days (model II), where /£ s = (1 + «d) ■ 
[D^D S /(D <j s Gpc)] with Dd, Dds and D s being the angular 
diameter distances between observer- lens, lens-source and 
observer-source, respectively and Zd the redshift of the lens. 
For typical lens (0.5) and source (1.5) redshifts, the delay is 
around 7//iso days (flat universe with O m = 1). The velocity 
dispersion of Gl is 97.6^5.4 • /d, s km s _1 (model II), where 
/d s = D s / Dds- The errors indicate the 99 per cent con- 
fidence interval, inferred from the Monte-Carlo simulations. 
The 68 per cent (la) confidence intervals are ~2.5 times 
smaller. Both the velocity dispersion and time delay depend 
on the chosen cosmological model through the angular di- 
ameter distances. A good description of the dependence of 
the model time delay and the angular diameter distances on 
the cosmologcal model is given in Helbig (1997). In a flat 
universe (fi m + Qa = 1) the difference in time delay is ^10 
per cent between O m = 1 and f2 m = 0.3. For a non-flat 
universe significant (^>10 per cent) differences in time delay 
are possible, depending on the dumpiness of matter and the 
combination of Qa (normalised cosmological constant) and 
f2 m (matter density). All other parameters in Table 4 are 
dimensionless and independent of the cosmological model. 

We conclude that model II is in substantial agree- 
ment with all available radio and optical observations of 
B1127+385. However, more detailed studies of the SB distri- 
butions of Gl and G2 are necessary to improve the models 
and tighten the confidence intervals. Moreover both lens and 
source redshifts are unknown. 



4.5 Galaxy colours, luminosities and M/L ratios 

Having convinced ourselves that Gl and G2 are the lens 
galaxies, we compare their colours (F555W-F814W«1.9- 
2.0) with the synthesized galaxy colours in Fukugita, Shi- 
masaku & Ichikawa (1995). For the different galaxy types 
approximate photometric redshifts of 0.3 (E), 0.4 (SO), 0.4 
(Sab), 0.7 (Sbc) and 0.9 (Scd) arc found. Gl and G2 could 
therefore be early type galaxies (including Sab) at low red- 
shift (0.3-0.4) or late type galaxies at high redshift (0.7-0.9). 
The integrated luminosities of Gl and G2 in B band for the 
different types of galaxies are: log 10 (Z/B©) ~ 8.4 — 2 log(/iso) 
(E), 8.8-21og(h 5 o) (SO and Sab), 9.6-21og(ft 5 o) (Sbc) and 
9.8-21og(/i 50 ) (Scd) for Gl (H = 50 • hso km s" 1 Mpc" 1 ), 
where we used the photometric redshifts found above and 
the B-F814W colours and K corrections from Fukugita et 
al. (1995). For G2 these values are 0.4 lower. 

At a redshift Zd = 0.3, 7=22.5 (Gl) corresponds to an 



Table 4. The minimum— \ 2 mass model parameters for 
B1127+385. Listed are the SD axis ratio (b/a)s, the velocity dis- 
persion an , the position angle (PA) of the major axis and the 
centre (x,y)Gl of Gl and the position and velocity dispersion 
of G2. The errors on axis ratio, velocity dispersion and position 
angle indicates the 99 per cent confidence interval inferred from 
Monte-Carlo simulations (Fig. 6). (x, 2/) src i/2 give the source po- 
sitions of the subcomponents A1&B1 and A2&B2, respectively. 
Furthermore, r^/ B and /iA,B give the derived flux density ratios 
and image magnifications. The coordinate system has been cen- 
tered on image Al. The definitions of and /? are given in 
Section 4.4.2. 





Model I 


Model II 




0.56 


r, rrA+0. 10 
°- 74 -0.12 


(x,S/)gi ( mas ) 


(-228.0, 


(-228.0, 




-38.0) 


-38.0) 


af 1 (kms" 1 ) 


108.9 ■ /| s 


"•<••;:! ■ / d a . 


PA G i (°) 


-26.4 


4 0+28.2 
^■ y -22.4 


(x,y)G2 (mas) 




(-749.0, 






-337.0) 


af 2 (kms- 1 ) 




0.79 ■ af 1 


(x,j/) 8 rcl (mas) 


-350.20±0.03, 


-439.81±0.03, 




-1.10±0.07 


-113.94±0.02 


(z,j/) src 2 (mas) 


-350.38±0.03, 


-439.89+.0.03, 




-1.62±0.07 


-114.42±0.02 


r Bl/Al 


0.56 


0.71 


r B2/A2 


0.57 


0.71 


MA1.A2 


-3.73 


-4.65 




-3.66 


-4.60 


MB1.B2 


2.09 


3.28 




2.09 


3.28 


At(Bi-Ai) (days) 


2.64 ■ 

J d,s 


1 ci +0.65 rb 

1 - Oi -0.60 Idy 


X 2 /NDF 


4.5 


1.1 



absolute magnitude Mj ~ —19 — 51og(ft 5 o). For an E and 
SO type galaxy this would mean it is ~4 mag underlumi- 
nous compared to E and SO type galaxies in the Hubble 
Deep Field (Mobasher et al. 1996). Placing the galaxy at 
higher redshifts would make the V — I colours of Gl and G2 
inconsistent with those of E or SO type galaxies (Fukugita 
et al. 1995). The absolute / magnitude is consistent however 
with somewhat later type spiral galaxies at higher redshifts 
(zd 0.5). This would also explain why the velocity disper- 
sions of Gl and G2 appears significantly smaller than those 
expected for L* E and SO type galaxies (e.g. Kochanek 1993, 
1994). 

Using the velocity dispersions listed in Table 4 (model 
II), the mass-to- light ratios (using the mass inside the Ein- 
stein radius, the photometric lens redshifts and the total B 
luminosity) of Gl and G2 are: (M/L B ) ~ 60 • ft 50 M /L b ,q 
(E), - 35 • h so (SO and Sab), - 10 • h 50 (Sbc) and ~ 10 • h 50 
(Scd) (assuming z s = 1.5). All of these mass-to-light ratios 
are significantly larger than normal galaxies of similar type, 
except for the higher redshift late-type spiral galaxies. Dust 
obscuration could increase the mass-to-light ratio, however 
it would also make the F555W-F814W colours less reason- 
able. 

We should note however that at intermediate redshift 
an error of 0.2 in redshift introduces ~1 mag error in B and 
as mentioned, dust obscuration (e.g. B1600+434; Jaunsen 
& Hjorth 1997; Koopmans, de Bruyn & Jackson 1998) or 
luminosity evolution (e.g. Bender, Ziegler & Bruzual 1996; 
Hudson et al. 1998) have not been taken into account. 
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Figure 6. From upper left to lower right: Monte-Carlo probability density distributions (of models with \ 2 < 11.3; 99 per cent confidence 
interval) of c^Gli (V a )s,Gi ar >d P-A.gi- The shaded regions (99 per cent confidence interval) indicate these parameters as derived from 
the SB distribution of Gl. The lower right panel shows x 2 as function of (<X|| j G2/ '||,Gl)- The observationally derived ratio falls well 
within the region containing 90 per cent of the \ 2 distribution. 



The results above therefore give only an indication. 
Redshifts of Gl and G2, and more accurate colours are vital 
to distinguish between galaxy types and their mass-to-light 
ratios. To obtain the redshifts of these galaxies in a rea- 
sonable integration time, one requires an 8 or 10-m class 
telescope (e.g. the Very Large Telescope (VLT) or Keck). 



4.6 Source 

We estimate I ~24.5 for the optical emission associated with 
radio component A. Correcting for a magnification of 4.6 
(Table 4), this corresponds to an intrinsic I~26, a luminos- 
ity of log 10 (Lio) ~ 9.8 — 21og(/i5o) at z s = 1.5 and an ab- 
solute magnitude Mj ~ —19 — 51og(fo 5 o) at z B = 1.5 (no 
K-corrections applied). 



Comparing this with the absolute I magnitudes of 
galaxies in the Hubble Deep Field, the source is most con- 
sistent with a spiral galaxy of type Sbc or later, as most 
E and SO type galaxies have Mi Si —23 (Mobasher et 
al. 1996). At a redshift of ~3 the luminosity increases to 
log 10 (Li Q ) ~ 10.5 — 21og(/i5o), still in the range of spiral 
galaxies (type Sab or later; Mobasher et al. 1996). We should 
note here that these values have not been corrected for evo- 
lution or dust absorption, internally or by the lens galaxies 
(Gl and G2). 

However, the CLASS lens systems B0712+472 (Jackson 
et al. 1998) and B1933+503 (Sykes et al. 1998; Jackson, pri- 
vate communication) also appear to have very low luminos- 
ity sources, which could indicate that a significant fraction 
of the weak (few mjy) flat spectrum radio sources are asso- 
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Figure 5. The 99 per cent confidence contour of the central SD 
position of Gl, determined only on the basis of the VLB A radio 
images. The cross gives the 99 per cent confidence region of the 
SB distribution of Gl relative to the optical component west of it. 
The two circles indicate the two radio components, assuming they 
are associated with the optical emission features (see text). With 
this assumption, the SB centre of Gl is in excellent agreement 
with its inferred SD centre. 

ciated with low luminosity objects, possibly late-type spiral 
galaxies. 



interval of the probability density distributions found from 
Monte-Carlo simulations. 

Having constructed a consistent model within the lens- 
ing hypothesis - explaining both the available radio and op- 
tical data of B1127+385 - it appears that B1127+385 most 
likely is a new gravitational lens system. 
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5 CONCLUSIONS 

A new gravitational lens system with two images separated 
by 701 ± 1 mas has been discovered in the CLASS survey. 
The two radio components have a flat spectrum between 
0.327 GHz (WSRT; WENSS), 1.7 GHz (MERLIN), 5.0 GHz 
(MERLIN) and 8.4 GHz (VLA). VLB A observations show 
substructure in both images. An HST /-band image reveals 
two emission features close to the radio components, which 
we identify with two lens galaxies. The colours and mass- 
to-light ratios of these galaxies seem to favour two late-type 
spiral galaxies at relatively high redshifts [zd 0.5). 

The VLBA radio structure and optical HST emission 
are consistent with a two-lens mass model, where a SIE and 
SIS mass distribution are placed on the SB centres of Gl 
and G2, respectively. This model is able to reproduce the 
separations and position angles of the VLBA substructure 
in radio components A and B, and their flux density ratios. 
Our best model has a reduced x 2 of 1.1. Assuming both 
lens galaxies are spiral galaxies, \ 2 slightly increases to 1.6. 
Omitting G2 in the mass model increases the reduced x 2 
to 4.5, which can therefore be excluded as an appropriate 
model with 99.6 per cent confidence. 

Our best model gives a SD axis ratio of 0.74^° '^i posi- 
tion angle of (4. 91^2 '4)° anu velocity dispersion of 97.6I5 4 • 
jj s km s -1 . The predicted time delay between radio compo- 
nents A and B is 1.5llQ go ■ fd B days. For a typical lens (0.5) 
and source redshift (1.5), a time-delay of ~7 days is expected 
for Ho=50 km s _1 Mpc -1 and f2 m = 1 in a flat universe. 
WSRT monitoring data is in hand, to see if the lensed object 
in B1127+385 is variable and therefore useful for determin- 
ing a time-delay between the lensed images. This time-delay 
can constrain the Hubble parameter (Refsdal 1964). The er- 
rors on these parameters indicate the 99 per cent confidence 
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